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Abstract 

This paper contains a phenomenological description of the whole 
U.S. forward rate curve (frc), based on an data in the period 1990- 
1996. We find that the average frc (measured from the spot rate) 
grows as the square-root of the maturity, with a prefactor which is 
comparable to the spot rate volatility. This suggests that forward 
rate market prices include a risk premium, comparable to the prob- 
able changes of the spot rate between now and maturity, which can 
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be understood as a 'Value-at-Risk' type of pricing. The instantaneous 
FRC however departs form a simple square-root law. The distortion 
is maximum around one year, and reflects the market anticipation of 
a local trend on the spot rate. This anticipated trend is shown to 
be calibrated on the past behaviour of the spot itself. We show that 
this is consistent with the volatility 'hump' around one year found by 
several authors (and which we confirm) . Finally, the number of inde- 
pendent components needed to interpret most of the FRC fluctuations 
is found to be small. We rationalize this by showing that the dynam- 
ical evolution of the FRC contains a stabilizing second derivative (line 
tension) term, which tends to suppress short scale distortions of the 
PRC. This shape dependent term could lead, in principle, to arbitrage. 
However, this arbitrage cannot be implemented in practice because of 
transaction costs. We suggest that the presence of transaction costs 
(or other market 'imperfections') is crucial for model building, for a 
much wider class of models becomes eligible to represent reality. 
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1 Introduction 



Finding an adequate statistical description of the dynamics of financial assets 
has a long history which is probably reaching a climax, now that enormous 
sets of time series are readily available for analysis. The case of the interest 
rate curve is particularly complex and interesting, since it is not the random 
motion of a point, but rather the history of a whole curve (corresponding to 
different maturities) which is at stake. The need of a consistent description 
of the whole interest rate curve is furthermore enhanced by the rapid devel- 
opment of interest rate derivatives (options, swaps, options on swaps, etc) 
[|I|]. Present models of the interest rate curve fall into two categories 0: the 
first one is the Vasicek model and its variants, which focuses on the dynamics 
of the short term interest rate, from which the whole curve is reconstructed 
(see for a recent review). The second one, initiated by Heath, Jarrow 
and Morton @] (see also [^), takes the full forward rate curve as dynamical 
variables, which is driven by one (or several) continuous Brownian motion, 
multiplied by a maturity dependent scale factor. Most models are primarily 
motivated by their mathematical tractability rather than by their ability to 
describe the data. For example, the fluctuations are often assumed to be 
gaussian, thereby neglecting 'fat tail' effects. More importantly, as we shall 
discuss below, the empirical shape of the interest rate curve can only be cap- 
tured by these models if one includes an additional parameter, the so called 
'market price of risk', which takes a rather large value. 

The aim of this paper is threefold. We first present an empirical study 
of the forward rate curve (frc), where we isolate several important features 
which a good model should be asked to reproduce. Rather than aiming 
at a precise statistical determination of parameters (which is known to be 
particularly difficult for all financial data), we wish to offer a qualitative 
description of the dynamics of the frc and try to decipher the information 
contained in its shape, in terms of an 'anticipated risk' and an 'anticipated 
bias'. Based on empirical evidence, we argue that the market fixes the interest 
rate curve through a Value-at-Risk type of condition, rather than through 
an averaging procedure, which is the starting point of the classical models 
alluded to above. Furthermore the 'anticipated bias' is found to be strongly 
correlated with the past trend on the spot rate itself. We then present a 
general class of string models, inspired from statistical physics, which describe 
the motion of an elastic curve driven by noise, and discuss how these models 
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offer a natural framework to account for some of the empirical results, in 
particular the small number of independent factors needed to describe the 
evolution of the FRC. Finally, we discuss the general concept of arbitrage and 
its relation with model building. We suggest that the argument of absence of 
arbitrage opportunities, usually used to restrict possible models of reality, is 
much too strong: in the presence of transaction costs and/or residual risks, 
qualitatively new classes of models can be considered. 

As an important preliminary remark, we want to stress that the adequacy 
of a model can be assessed from two rather different standpoints. The first 
one, which we shall adopt here, is to see how well a given model describes the 
dynamics of the primary object, namely the forward rate curve. The second 
one, more concerned with derivative markets, asks how consistent are the 
determination of the model's parameters across several different derivative 
products, irrespective of the ability of these parameters to reproduce the 
statistics of the underlying asset. A well known example is the case of the 
Black-Scholes theory of options: while there should be no volatility smile in 
a Black-Scholes world, it is still an intersting question from an applied point 
of view to know whether the implied volatility on plain vanilla options can 
be used to price (within a Black-Scholes scheme) exotic options. We will 
however not consider in the present paper the consequence of our model to 
the problem of derivative pricing. 

2 Statistical analysis of the forward rate curve 

2.1 Presentation of the data and notations 

The forward interest rate curve (frc) at time t is fully specified by the 
collection of all forward rates f{t,9), for different maturities 9. It allows 
for example to calculate the price B{t, 9) at time t of a (zero-coupon) bond, 
which pays 1 at time t -\- 9. It is given by: 



r(t) = f{t.9 = 0) is called the 'spot rate'. Note that in the following 9 is 
always a time difference; the maturity date T is t + 9. 




(1) 
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2.2 The data set 



Our study is based on a data set of daily prices of Eurodollar futures con- 
tracts. The Eurodollar futures contract is a futures contract on an interest 
rate, as opposed to the Treasury bill futures contract which is a futures 
contract on the price of a Treasury bill. The interest rate underlying the 
Eurodollar futures contract is a 90-day rate, earned on dollars deposited in 
a bank outside the U.S. by another bank. When interest rates are fixed, a 
well known arbitrage argument implies that forward and futures contracts 
must have the same value. But when interest rates are stochastic in principle 
forward contracts and future contracts are no longer identical - they have dif- 
ferent margin requirements - and one may expect slight differences between 
futures and forward contracts , which we shall neglect in the following. We 
have furthermore checked that the effects discussed below are qualitatively 
the same as those observed when FRC is reconstructed from swap rates. The 
interest of studying forward rates rather than yield curves is that one has a 
direct access to a 'derivative' (in the mathematical sense - see Eq.(||)), which 
obviously contains more precise informations. 

In practice, the futures markets price three months forward rates for 
fixed expiration dates, separated by three month intervals. Identifying three 
months futures rates to instantaneous forward rates, we have available time 
series on forward rates /(t,Tj — t), where Tj are fixed dates (March, June, 
September and December of each year), which we have converted into fixed 
maturity (multiple of three months) forward rates by a simple linear interpo- 
lation between the two nearest points such that Ti —t < 9 < Tj+i — t. Days 
where one contract disappears are not included, to remove possible artefacts. 
Between 1990 and 1996, we have at least 15 different Eurodollar maturities 
for each market date; longer maturities are in general less liquid than short 
maturities; we however believe that such a difference in liquidity does not af- 
fect the qualitative conclusions drawn in this paper (this is partly confirmed 
by our analysis of the swap rates). Between 1994 and 1996, the number of 
available maturities rises to 30 (as time grows, longer and longer maturity 
forward rates are being traded on future markets); we shall thus often use 
this restricted data set. Since we only have daily data, our reference time 
scale will be r = 1 day. The variation of /(t, 9) between t and t + t will be 
denoted as 6f{t,9): 

6f{t,9) = fit + r,9)-fit,9) (2) 
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It would obviously be very interesting to extend the present analysis to intra- 
day fluctuations. 

2.3 Quantities of interest and data analysis 

The description of the FRC has two, possibly interrelated, aspects : 

(i) what is, at a given instant of time, the shape of the FRC as a function 
of the maturity 6 ? 

(ii) what are the statistical properties of the increments 5f{t, 6) between 
time t and time t + r, and how are they correlated with the shape of the PRC 
at time t ? 

The two basic quantities describing the FRC at time t is the value of 
the short term interest rate f(t, 6'min) (where 6'min is the shortest available 
maturity), and that of the short term/long term spread s(t) = /(t,6'max) — 
f{t, dmin), where 6'max is the longest available maturity. The two quantities [] 
r(t) ~ fit, 9mm), s(t) are plotted versus time in Fig. 1; note that 

- The volatility a of the spot rate r{t) is equal to O.S%/ y/jeai. This is 
obtained by averaging over the whole period. However, as previously noticed 
[§, P], there seems to be a systematic correlation between a and r, which can 
be parametrized as o" oc r^. For the period considered (1990-96), we find a 
much smaller value /? ~ 0.4 than the one reported in |^ (3 = 1.5, based one 
monthly data from 1964 to 1989, or /3 = 1.1 based on weekly data from 
1984 to 1993. 

- s{t) has varied between 0.53% and 4.34%. Contrarily to some european 
interest rates on the same period, s{t) has always remained positive, a sit- 
uation we shall refer to as 'normal'. (This however does not mean that the 
FRC is increasing monotonously, see below). 

In order to rationalize the empirical data, we shall postulate that the 

^We shall from now on take the three month rate as an approximation to the spot rate 
rit). 

^The dimension of r should really be % per year, but we conform to the habit of 
quoting r simply in %. Note that this can sometimes be confusing when checking the 
correct dimensions of a formula. 
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Figure 1 
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Figure 1: The historical time series of the spot rate r(t) from 1990 to 1996 
(top curve) - actually corresponding to a three month future rate (dark line) 
and of the 'spread' s(t) (bottom curve), defined with the longest maturity 
available over the whole period 1990-1996 on future markets, i.e. Omax = four 
years. (In restricted period (94-96), this maturity grows to 6'max = 8 years.) 
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whole FRC oscillates around an average shape, and parametrize f{t, 9) as 

fM=r{t)+s{t)^^ + am (3) 

where F is a certain (time independent) function and .^(t, 9) represents the 
deviation from the (empirical) average curve, in the sense that 

mo)) = ^ (4) 

where (...) represents a time average. By definition, ^(t, 6*111111) = ^(^j^'max) = 
0. 

Fig. 2 shows the scaling function Y{u), determined by averaging the 
difference f{t,9) — r{t) over time. Interestingly, the function Y{u) is rather 
well fitted by a simple ^/u law. This means that on average, the difference 
between the forward rate with maturity 9 and the spot rate is equal to 
A\/9, with a proportionality constant A = 0.85%/iyyear which turns out 
to be nearly identical to the spot rate volatility. We shall propose a simple 
interpretation of this fact below. 

Let us now turn to an analysis of the fluctuations around the average 
shape. Fig. 3 shows for different instants of time the whole 'deviation' ^{t, 9) 
as a function of 9. 

These fluctuations are actually similar to that of a vibrated elastic string. 
The average deviation A{9) can be defined as: 

A{9) := ^mer) (5) 

which is plotted in Fig. 4, for the period 94-96. The maximum of A is 
reached for a maturity of 6** =1 year. 

A more precise study of the fluctuation modes consists in studying the 
equal time correlation matrix A4{9,9') defined as: 

Mi9,9') = {at,9)at,9')) (6) 

Note that M{9,9) = A{9)'^. The eigenvalues of this matrix are plotted 
in Fig. 5, as a function of their rank q = 1, 2, in log- log coordinates. We 



■^The relation between this decomposition and the more standard Principal Component 
Analysis is clarified below. Note in particular that r{t), s{t) and ^{t,9) are a priori not 
independent. 
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Figure 2: The average PRC in the period 94-96, as a function of the maturity 
Q. We have shown for comparison a one parameter fit with a square-root law, 
A{\fQ — V^min)- The same behaviour actually extends up to ^max = 10 
years, which is available in the second half of the time period. 
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Figure 3: This shows the deviation ^{t, 6) as a function of 6 for different times 
t. Note that this resembles the motion of an elastic string, with a maximum 
around 6 — 1 year. 
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Figure 4: Root mean square deviation A(^) from the average FRC as a func- 
tion of 6. Note the maximum for 9* = 1 year, for which A ~ 0.38%. 
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have shown for comparison a and decay (see section 5), indicating 
that these eigenvalues are decreasing very fast with q. We have furthermore 
calculated the same eigenvalues in the presence of some 'artificial' noise (i.e. 
a random variable of zero mean and width 0.04% added to each /(t, 9) inde- 
pendently). One sees that it only affects the high q modes, while leaving the 
modes with g < 9 relatively stable, and which are thus statistically meaning- 
ful. The deviation from the average shape can be written as: 

at.e) = Y.M,ut)'^,{e) (7) 

where {C,q{t)^q'{t)) = Sq^q>. Exactly as an elastic string, the first mode has 
no nodes (and is found to be nearly proportionnal to A(0)), while the second 
\l/2 has one node. 

2.4 Statistics of the daily increments 

We now turn to the statistics of the daily increments Sf{t, 9) of the forward 



rates, by calculating their volatility a{6) = y {6f{t, 0^), their excess kurtosis 

<0) - ^4(^) - 3 (8) 
and the following 'spread' correlation function: 

= ^^Ho-) 

which measures the influence of the short term interest fluctuations on the 
other modes of motion of the FRC. 

Fig. 6 shows a {6) and k{6). Somewhat surprisingly, cr(0), much like A(^) 
has a maximum around 6* = 1 year, a feature already noticed in, e.g., []TU|, |TT|]. 
The order of magnitude of a{9) is 0.05%/A/day, or 0.8%/^year. The kurtosis 
k{6) is rather high (of the order of 5), and only weakly decreasing with 6. 

Finally, C{6) is shown in Fig. 7; its shape is again very similar to those of 
A(^) and cr(^^), with a pronounced maximum around 6 = 1 year. This means 
that the fluctuations of the short term rate are amplifiedior maturities around 
one year. We shall come back to this important point below. Note that 
C{6) goes to zero for large maturities, which means that the short term rate 
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Figure 5: The eigenvalues of the instantaneous fluctuation matrix 
Ai{6,6'), as a function of their rank q = 1,2, in log-log coordinates. We 
have shown for comparison a and a decay (see section 5), indicat- 
ing that these eigenvalues are decreasing very fast with q. Also shown for 
comparison the eigenvalues of Ai {6, 9') when an artificial noise (uniform in 
[—0.02, -1-0.02]) is independently added to each point of the FRC, to estimate 
the relevance of the information contained in the high q eigenvalues. 
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Figure 6: The daily volatility and kurtosis as a function of maturity. Note 
the maximum of the volatility for 6 = 6*, while the kurtosis is rather high, 
and only very slowly decreasing with 6. The two curves correspond to the 
periods 90-96 and 94-96, the latter period extending to longer maturities. 
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and the short term/long term spread evolve, to a first approximation, in an 
uncorrelated manner. A simple fit of C{6), shown in Fig. 7, is given by C{6) = 
c6exp{—r9), with c = 1.5 year~^ and F = 1.0 year^^. An interpretation of 
the values of c and F will be given below. 

Much as above, one can define a correlation matrix for the variations of 
^{t,e) as: 

^^i9,9') = {5at,e)5at,9')) (lo) 

Its diagonalisation gives results very similar to the diagonalisation of A^, in 
particular concerning the fall-off of the eigenvalues as a function of q and the 
shape of the low q modes. 

Owing to the rapid decrease of the eigenvectors of M. and A^, a good 
approximation thus consists in retaining only the so called 'butterfly' mode 
of deformation and to write 0: 



f(t,9) = r{t) + s{t)y 



9 



9r. 



(11) 



where y{u) := a/u- That only a rather small degrees of freedom are needed 
to describe most of the FRC's fluctuations was already discussed on several 
occasions [O, |T( 



13, 14], although not exactly in the present terms. 



3 Classical models 



3.1 Vasicek 



The simplest FRC model is a one factor model due to Vasicek ||T^, where the 
whole term structure can be ascribed to the short term interest rate, which 
is assumed to follow a stochastic evolution described as: 

dr{t) = n{ro - r{t)) + adW{t) (12) 

where tq is an 'equilibrium' reference rate, Q describes the strength of the 
reversion towards tq (and is the inverse of the mean reversion time), and 
dW{t) is a Brownian noise, of volatility 1. In its simplest version, the Vasicek 
model prices a bond maturing at T as the following average: 

B{t, T) = ( exp - y du r{u) ) (13) 
"^We have included the factor Mi in the definition of ^I^i 
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Fig 7 
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Figure 7: The correlation function C{9) (defined by Eq. ||) between the daily 
variation of the spot rate and that of the forward rate at maturity 6, in the 
period 94-96. Again, C{6) is maximum for 6 = 6*, and decays rapidly beyond. 
A simple fit using the function c6 exp{—r6) is shown for comparison, which 
leads to c = 1.5 year~^ and T = 1.0 year~^. Note also that C(6'max) is close 
to zero, indicating that, in a first approximation, the evolution of the spot 
rate and of the spread are independent. 
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where the averaging is over the possible histories of the spot rate between 
now and the maturity, where the uncertainty is modelled by the noise W. 
The computation of the above average is straightforward and leads to (using 
Eq. (I): 

fit, 6) = r{t) + (ro - r(t))(l - e"™) - - e-^^ (14) 

The basic results of this model are as follows: 

• Since (ro — r(t)) = 0, the average of f{t, 6) — r{t) is given by 

(/(t, 6) - r{t)) = -a'/2n\l - e-™)^ (15) 

and should thus be negative, at variance with empirical data. Note that in 
the limit fi^ <^ 1, the order of magnitude of this (negative) term is very 
small: taking a = 1%/^year and 9 = 1 year, it is found to be equal to 
0.005% ! 

• The volatility a{6) is monotonously decreasing as exp —Q6, while the 
kurtosis k{9) is identically zero (because W is gaussian) 

• The correlation function C{6) is proportionnal to the volatility cr(^), 
and thus does not exhibit the marked maximum shown in Fig. 7. 

• The variation of the spread s(t) and of the spot rate should be perfectly 
correlated, which is not the case (see Fig. 7): more than one factor is in any 
case needed to account for the deformation of the FRC. 

While it is easy to make k non zero by taking a discrete time process 
where dW is non gaussian, it is awkward to account for a maximum in the 
volatility (j{6) within such an approach (see the discussion below). However, 
the most obvious inconsistency of this model is the fact that the average 
spread should be negative. A way out is to introduce the 'market price 



of risk': as shown by Vasicek |jT5|, the probability measure over which the 
average (0) is performed is not necessarily the historical average. Arbitrage 
arguments allow a 'change of measure', which in the present case simply 
amounts to correcting the 'true' (i.e historical) value of tq to an effective 
value ro + Act, where A is the market price of risk. With this correction, one 
finds that: 

(/(t,^)-r(t)) = A(T(l-e-™) (16) 

A fit of the empirical data with such a formula leads to = 6.17 10~^/year 
and Xa = 2.26%. The value of Q corresponds to a mean reverting time of 
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Figure 8: Fit of the FRC with a Vasicek term structure, Eq. (|1^) with a non 
zero market price of risk. This is a two parameter fit - compare to the one 
parameter fit in Fig. 2 - which is however found to be acceptable. One finds 
Vt = 6.17 10^^/year and Xa = 2.26%. Note that we neglect the difference 
between r{t) and /(t,6'min). 
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around 16 years, which is worrying since the data set is only 7 years long - 
it is thus not really consistent to set {tq — r{t)) to zero. The market price 
of risk q corresponds to demanding for a one year maturity bond an extra 
return (over the spot rate) of around 0.15%. 



3.2 Hull and White 

An interesting extension of Vasicek's model designed to fit exactly the 'initial' 
FRC f(t = 0, 9) was proposed by Hull and White [0. It amounts to replacing 
the above constants Q and ro by time dependent fonctions. For example, 
ro(t) represents the anticipated evolution of the 'reference' short term rate 
itself with time. These functions can be ajusted to fit f{t = 0,6) exactly. 
Interestingly, one can then derive the following relation (for a zero market 
price of risk A): 

(^)h|m,) (17, 

up to a term of order which turns out to be negligible, exactly for the 
same reason as explained above. On average, the second term (estimated by 
taking a finite difference estimate of the partial derivative using the first two 
points of the FRc) is definitely found to be positive, and equal to 0.8 %/year. 
On the same period (90-96), however, the spot rate has decreased from 8.1 
% to 5.9 %, instead of growing by 5.6%. 

In simple terms, both the Vasicek and the Hull- White model mean the 
following: the FRC should basically reflect the market's expectation of the 
average evolution of the spot rate (up to a correction on the order of a^, 
but which turns out to be very small - see above). However, since the 
FRC is on average increasing with the maturity (situations when the FRC is 
'inverted' are comparatively much rarer), this would mean that the market 
systematically expects the spot rate to rise, which it does not. It is hard 
to believe that the market persists in error for such a long time. Hence, 
the upward slope of the FRC is not only related to what the market expects 
on average, but that a systematic risk premium is needed to account for this 
increase. Within a classical framework, this is attributed to the 'market price 
of risk', which is however not a directly measurable quantity. In section (^), 
we will propose a more direct interpretation of this risk premium. 

Let us however mention that some empirical results can be accomodated 
by the Hull- White model and its extensions |T^, 0, |1^. For example, a 
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'humped volatility' can be obtained by choosing f2 to be time dependent in 
such a way that it starts being negative at time t = (meaning that somehow 
the spot rate should escape from its reference value tq) until a certain time 
To beyond which it remains positive. It is easy to show in that case that 
the volatility is maximum for a maturity 6 = Tq. However, apart from the 
fact that the interpretation of this curious shape for Q{t) is not clear, one 
can also show that in that case the average FRC should display an inflexion 
point at ^ = To, which is not seen empirically. In order to fit independently 
the shape of the FRC and the shape of the volatility, one should turn to two- 
factor Hull- White models [0, in which each one of the two mean reversal 
process can be more easily interpreted (see our own discussion below). 



3.3 Heath- Jarrow-Morton 

A second, more recent, line of thought, consists in writing stochastic dif- 
ferential equations for each of the forward rates. In this case, today's PRC 
is by construction described exactly, and provides the initial condition for 
the stochastic evolution. The simplest, time translation invariant, one factor 
model reads P, H : 



df{t,9) 



df{t,9) , ■ 



dt + a{e)dW{t) (18) 



with a certain maturity dependent volatility o"(6'). If the expiration date 
T = t + 6 rather than the maturity is kept fixed, the corresponding rate 
/(t, T) = /(t, ^ = T - t) is such that: 

df{t, T) = df{t, 6) - ^^^^dt = n{e)dt + a{e)dW{t) (19) 

Ou 

Absence of arbitrage opportunities then impose, within this continuous time 
framework, a relation between /i and a which reads: 

^(9) = a{d) f a{e')de' + Act(^), (20) 
Jo 

where A is again the market price of risk. This model has primarily been 
devised to price consistently (within a no-arbitrage framework) interest rates 
derivatives rather than to represent faithfully the historical evolution of the 
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FRC itself. It is however interesting to remark that the above criticism still 
holds if one wants to interpret data using Eq. (^): since the average value of 
the slope of the FRC can only be generated by the maturity dependent drift 
IJ,{9), a non zero value of A is needed in order to reproduce its correct order 
of magnitude - the order cr^ term is again much too small, and can for all 
practical purposes be set to zero (this was first noticed in pO[). Furthermore, 
Eq. (^) does not ensure that the shape of the PRC remains 'realistic' with 
time: one the contrary, the FRC tends to distort more and more as time 
passes when the volatility has a non trivial maturity dependence. We shall 
discuss in section 5 below how the introduction of a 'line tension' may heal 
this disease. 

In conclusion, the classical models are only consistent with empirical data 
provided a rather large market price of risk is included. In technical terms, 
this means that, in the context of arbitrage theories, the 'risk-neutral' prob- 
ability measure to be used for - say - derivative pricing cannot be identified 
with the empirical probability. Other theoretical difficulties inherent to ar- 
bitrage theories of the interest rate curve have been discussed in |p.8|| . 

4 Risk-premium and the \/9 law 

4.1 The average FRC and Value-at-Risk pricing 

The observation that on average the PRC follows a simple \/9 law (i.e. 
(/(t, 6*) — r{t)) oc \/9) suggests a more intuitive, direct interpretation. At 
any time t, the market anticipates either a future rise, or a decrease of the 
spot rate. However, the average anticipated trend is, on the long run, zero, 
since the spot rate has bounded fluctuations. Hence, the average market's 
expectation is that the future spot rate r{t) will be close to its present value 
r{t = 0). In this sense, the average PRC should thus be flat (again, up to a 
small cr^ correction). However, even in the absence of any trend in the spot 
rate, its probable change between now and t = 6 is (assuming the simplest 
random walk behaviour) of the order of cr\/^, where a is the volatility of the 
spot rate. Hence, money lenders are tempted to protect themselves against 
this potential rise by adding to their estimate of the average future rate a 
risk premium of the order of a\/9 to set the forward rate at a satisfactory 
value. In other words, money lenders take a bet on the future value of the 
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spot rate and want to be sure not to lose their bet more frequently than - 
say - once out of five. Thus their price for the forward rate is such that the 
probability that the spot rate at time t + 9, r(t + 9) actually exceeds f{t, 9) 
is equal to a certain number p: 

POO 

/ dr'P{r',t + 9\r,t) =p (21) 
Jf{t,e) 

where P{r' ,t'\r,t) is the probability that the spot rate is equal to r' at time 
t' knowing that it is r now (at time t). Assuming that r' follows a simple 
random walk centered around r{t) then leads to Q 

f{t, 9) = r{t) + AaVd A = V2 erfc"^(2p) (22) 

which indeed matches the empirical data, with p ~ 0.16. 

At this point, we thus depart both from the so-called 'unbiased expec- 
tation hypothesis' 0] or from Vasicek type of models, where it is assumed 
that today's FRC tells us something about the average future evolution of 
the spot rate, possibly corrected by an (unmeasurable) market price of risk. 
This, for example, is used to calibrate the time scale fi^^ which appears in 
the Vasicek model, Eq. (|T^. In our mind, the shape of today's FRC must 
rather be thought of as an envelope for the probable future evolutions of the 
spot rate. The market seems to price future rates through a Value at Risk 
procedure (Eqs. ^ |2^) rather than through an averaging procedure. In a 
strict sense, Eq. ( ^2|) is not acceptable from the point of view of arbitrage, 
since buying the spot rate and selling the forward rate would lead to a profit. 
However, one should keep in mind that this profit is not certain because Eq. 
( |22D only describes the PRC in an average sense; one should also take into 
account the random 'deformation' of the curve described by C,{t,9) (which 
we interpret below) which makes the above strategy risky. Nevertheless, the 
simple strategy of lending long-term money and borrowing short-term money 
is obviously a major source of income for most banks ! 



^This assumption is certainly inadequate for small times, where large kurtosis effects 
are present. However, one the scale of months, these non gaussian effects can be considered 
as small 
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4.2 The instantaneous PRC. The anticipated bias. 

Let us now discuss, along the same lines, the shape of the PRC at a given 
instant of time, which of course deviates from the average square root law. 
As discussed in section 2, the deviation ^(t, 6) is actually well approximated 
by the first 'butterfly' mode which has the shape drawn in Fig. 4. 

We interpret this as follows: for a given instant of time t, the market actually 
expects the spot rate to perform a biased random walk. We shall argue that 
the market estimates the trend m{t) by extrapolating the past behaviour of 
the spot rate itself. However, it appears that the market also 'knows' that 
trends on interest rates do not persist for ever in time; the 'anticipated' bias 
is thus in general maturity dependent. Hence, the probabihty distribution 
P{r',t + 6\r,t) used by the market is not centered around r{t) but rather 
around: 

rt+e 

r{t) + j dum{t,t + u) (23) 

where m{t, t') can be called the anticipated bias at time t', seen from time t. 

As discussed above, it is reasonable to think that the market estimates m 
by extrapolating the recent past to the nearby future. Mathematically, this 
reads: 

POO 

m(t,t + u) = mi(t)g(u) mAt):= dv K(v)5r(t - v) (24) 

JO 

where K{v) is an averaging kernel of the past variations of the spot rate. 
We will call g{u) the trend persistence function; it is normalized such that 
^(m = 0) = 1, and describes how the present trend is expected to persist in 
th efuture. Eq. (^) then gives: 

fit, 6) = r(t) + AaVe + mAt) f du g{u) (25) 

JO 

Identifying this expression with Eq. (p^), one finds that: 

• As discussed above, the spread s is associated to the volatihty of the 
spot rate. Note however that other sources of risk, such as the exchange rate 
fluctuations for foreign investors, might also be included in the 'effective' 
volatility of the spot rate used in Eq. (PT|). 

• The persistence function is simply the derivative of the first mode '^i{9), 
or A(6'). This function is plotted in Fig. 9; the interesting point - related to 
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the existence of a maximum in A{6) - is that it becomes negative (but small) 
beyond one year. This means that the market anticipates a trend reversion 
on the scale of one year. 

This interpretation furthermore allows one to understand why the three 
functions introduced above, namely A(^), a{6) and the correlation function 
C{6) have similar shapes. Indeed, taking for simplicity the kernel K{v) to be 
7exp[— 7v], one finds: 

dmi{t) = —'jmidt + '^dr{t) (26) 
Hence, the correlation function C{d) is nothing but: 

c{e) = -i f dug{u) (27) 

Jo 

(we have assumed in the above expression that C{9^^^ ~ 0, which is not a 
bad approximation - see Fig. 7). On the other hand, since from Eq. (p^) 
{m\) = 7CT^(0)/2 one also finds: 

m = i/|^(0) /J duGiu) = ^C{e) (28) 



thus showing that A(^) and C{6) are proportionnal. Actually, even the nu- 
merical prefactor predicted by this simplified description is quite good: using 
the value of 7 = c determined by the fit shown in Fig. 7, the predicted value 
of A(^*) is found to be ~ 0.29%, instead of the observed 0.38% (see Fig. 4). 
Turning now to the volatility cr(0), one finds that it is given by: 

a\e) = [1 + C{9)f a^O) + a^^y{9y (29) 

where is the contribution of the spread volatility. We thus see that the 
maximum of a{9) is indeed related to that of C{9). Intuitively, the reason for 
the volatility maximum is as follows: a variation in the spot rate changes that 
market anticipation for the trend mi(t). But this change of trend obviously 
has a larger effect when multiplied by a longer maturity. For maturities 
beyond one year, however, the decay of the persistence function comes into 
play and the volatility decreases again. 



Interestingly, the form (p9D suggests that the volatility should increase 



again for very large 9, due to the contribution of the last term. This can 
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Figure 9: The 'persistence' function Q{6) describing how the market antic- 
ipates the persistence in the future of the observed past trend of the spot 
rate. Note that Q{0) becomes negative beyond one year, meaning that the 
market anticipates a trend reversion on the scale of one year. 
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indeed be seen on a restricted data set which allows to reach ^ = 10 years 
(Fig. 10), where one sees that a{6) reaches a minimum around ^ = 8 years. 
This effect can also be related to the slow decay of the kurtosis with 9 shown 
in Fig. 6, which reflects the fact that short maturities are essentially sensitive 
to one factor (the spot rate fluctuations), while longer maturities are progres- 
sively affected by a second factor (the anticipated risk) - adding independent 
contributions indeed reduce the kurtosis. 

Taking for simplicity Jq du Q{u) = Oexp—TO, as suggested by Fig. 7, 
one can now fit the instantaneous FRC by the expression (|25|) , for example 
fixing F. One then extracts from the FRC the anticipated bias mi{t) as a 
function of t. We find that mi{t) is correlated with the past evolution of the 
spot rate, on the scale of one year (i.e. the past time scale determining mi, 
is comparable to the time scale of the persistence function, F^^). This 
is perfectly compatible with the value of 7 = c extracted directly from the 
fit shown in Fig. 7, which gives 7"^ = 8 months. 

Finally, we want to note that Eq. (126|) is actually very similar in spirit 



to the second equation in Hull and White's two factor model |T0[ , where the 
'noise' term in the second equation is actually the spot rate itself. This model 
was constructed in an ad-hoc way to reproduce the volatility hump and its 
intuitive meaning was not explicited. We believe that the above idea of an 
'autoregressive' anticipated biais is a rather natural financial interpretation. 



5 The forward rate curve as a vibrating string 

In the previous section, we have thus argued that the shape of the instanta- 
neous FRC is fixed by two distinct type of anticipations: 

- An anticipated bias, whose influence is expected to decay with time, 
and whose amplitude is determined by the recent evolution of the spot rate. 

- A anticipated risk, which gauges the potential motion of the spot rate, 
and which is added as a risk premium by money lenders. 

Hence, as noticed by several people, the whole FRC evolves according to 
rather few 'factors' while in principle each maturity 6 could feel the 

influence of some independent factor. However, one also expects that the 
'forces' determining the evolution of each points of the FRC act as to prevent 
the FRC from 'blowing apart' with time, and to diffuse the information from 
one maturity to the next. A natural mechanism is the following: the dynami- 
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Figure 10: Volatility cr as a function of 9, as in Fig. 6, but for a restricted 
time period allowing to extend the range of maturities to 10 years. A volatil- 
ity minimum appears around ^ = 8 years, beyond which the volatility of the 
spread becomes important. 
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cal evolution should be such that large deviations between nearby maturities 
are improbable - it is plausible that market forces tend to keep the forward 
rate f(t, T) close to its local average [f(t, T — e) + f(t, T + e)]/2, where here 
e = 3 months. The simplest choice is then to write: 

6 fit, T)=V [/>, T - e) + />, T + e) - 2/(t, T)] + r^{t, T) (30) 

where V measures the strength of this force ('line tension') which could in 
principle be 9 dependent. Eq. (|30D describes a mean reversion of the forward 
rate towards its local average. It is the dynamical equation of a string of 
beads connected by springs, subject to a random force term rj which can 
reflect, in particular, the variations of the forward rate which are specific to 
that given maturity (for example, the influence of a major political event 
which is scheduled around that particular date.). 

The term proportional to T) would become, in the continuum limit, a 
second order derivative d"^ f /dT'^, which is in principle not allowed from arbi- 
trage arguments, at least when the 'noise' ri{t,T) contains a finite number of 
independent components. Intuitively, this can be understood in the simple 
case where ri{t,T) is independent of T. Then, a possible riskless winning 
strategy is to buy (resp. sell) one contract of maturity T if 

fit, T - e) + fit, T + e) - 2/(t, T) > (31) 

(resp. < 0), and then to hedge the position by buying (or selling) the ap- 
propriate number of the shortest maturity contract. Then, because of the 
restoring 'string' force, the resulting gain will be positive. The result of such 
a strategy (used simultaneously on all maturities) is shown in Fig. 11, in the 
absence of transaction costs. This curve shows indeed a clear positive slope, 
thereby directly showing the existence of a V term affecting the temporal 
evolution of the FRC. The fact that the above strategy is not riskless is due 
to the presence of more than one random factor rj. More importantly, we 
have found empirically that the above growing P&L is completely destroyed 
by transaction costs. 

Therefore, the presence or absence of shape dependent terms (such as the 
above second order difference term) cannot be decided by arbitrage arguments 
in the presence of transaction costs. 

Still, the very existence of these terms is extremely important to under- 
stand the qualitative dynamical evolution of the FRC. In particular, the 'line 
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Figure 11: Empirical Profit and Loss curve corresponding to the simple 
'arbitrage' strategy proposed in the text, which exploits the presence of a 
'string' force proportionnal to the local second (discrete) derivative of the 
FRC. On average, this strategy clearly leads to a positive profit, although 
reasonable transaction costs ruins it completely. However, the presence of 
a non zero T) is very important since it allows one to understand why the 
effective number of factors influencing the FRC is small (see Eq. ( p4D below). 
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tension' ensures that FRC becomes (and remains) smooth even if the initial 
condition is not, of if the term structure of the volatihty is not. This is not 
the case within Heath- Jarrow-Morton hke formulations, where V = 0, and 
where temporal evolutions tend to be ill-behaved . 



We want to show here that the mere existence of a line tension allows one 
to understand qualitatively the fast decay of the eigenvalues of the correlation 
matrix, as observed in Fig. 5. If one introduces 4>it,T) = f{t,T) — f{t,t), 
and takes for notational simplicity the continuous time limit, one finds that 
evolves according to: 

T) = V^^^dt + rfr/(t, T) t) = (32) 

Defining the eigenmodes of the operator which vanish at T = t, i.e. : 
(f)q(t,T) = Mq{t) sin(g(T — t)), the evolution of Mq(t) is then given by: 

^ = -Vq'Mq + ^,{t) (33) 

where rjq is the Fourier transform of the noise term. In this form one sees 
that: 

• The 'wavevector' q is the index labelling the eigenmodes of the matrix 
M. introduced above. 

• The 'lifetime' of a perturbation of wavevector q is inversely propor- 
tionnal to Vq^. High g's (corresponding to short maturity differences) relax 
faster - this corresponds to the fact that if the FRC is distorted on very small 
scales, this will probably not persist in time. 

• The average (over time) amplitude in mode g, proportional to the g*^ 
eigenvalue of ^A is given by: 

= ^ (34) 

If Tjit, T) was independent for each maturity, (|?7gP) would be independent 
of g, already leading to a decay of (|Mgp), due to the restoring force V. 
The presence of correlations in the noise along the maturity axis acts to make 
this decay even faster. For example, if the correlations are exponentially 
decaying for different T's (i.e. as exp —K\T — T'|), one would find: 

« (35) 
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hence leading to a much faster still decay of (|Mgp). As shown in Fig. 5, the 
observed decay is indeed intermediate between and q~^. 

The main conclusion of this section is that a 'line tension' term indeed 
exists in the dynamical evolution of the FRC, acting as to reduce the distor- 
tions along the line. This line tension is responsible for reducing the effective 
number of factors needed to interpret the shape of the whole FRC, even when 
independent shocks affect each different maturities 0. Arbitrage arguments 
in the absence of transaction costs impose that V = 0, a. conclusion that 
transaction costs allows one to by-pass. 

6 Summary-Conclusion 

The main results contained in this paper, based on the study of the whole U.S. 
forward rate curve since 1990 (and which we confirmed on the corresponding 
swap rates), are the following: 

• The average PRC (measured from the spot rate) grows as the square- 
root of the maturity, with a prefactor which is comparable to the spot rate 
volatility. This strongly suggests that forward rate market prices include a 
risk premium, comparable to the probable changes of the spot rate between 
now and maturity. This interpretation seems to us more natural (although 
in the same spirit) than the one based on an unobservable 'market price of 
risk'. 

• The instantaneous FRC departs form a simple square-root law. The 
distortion is maximum around one year, and reflects the market anticipa- 
tion of the trend on the spot rate (as in more traditionnal models of interest 
rates pO|). Somewhat surprisingly, however, this trend appears to be cali- 
brated on the past behaviour of the spot itself: if the spot rate has decayed 
over the past months, the one year forward rate is below the average square 
root extrapolation, and vice-versa. This is consistent with the fact that the 
volatility is maximum for one year maturities. 

• The number of independent components needed to interpret most of the 
PRC fluctuations is rather small, although in principle a random contribution 
should be assigned to each maturity. We rationalize this finding by showing 

^The idea of projecting the FRC on a few smooth functions of the maturity can also be 
found in |2^, ahhough the basic mechanism underlying such a reduction, in terms of 
a line tension, was not discussed. 



31 



that the dynamical evolution of the PRC contains a stabilizing 'line tension' 
term, which tends to suppress short scale distortions of the frc. This shape 
dependent term in the dynamical evolution is not usually considered in in- 
terest rate models, because it leads, in principle, to arbitrage. However, this 
arbitrage cannot be implemented in practice because of residual risks and 
transaction costs. This is important because the long time statistical proper- 
ties of a line subject to random forcing is very different when the line tension 
term is present (even very small) or strictly zero: the 'line tension' term (as 
all second derivative terms) lead to a smoothing of the singularities, which 
otherwise survive or even develop in the course of the evolution. We thus 
conclude that the presence of transaction costs (or other market 'imperfec- 
tions', such as residual risk) is crucial for model building, for a much wider 
class of models becomes ehgible to represent reality. 

The present work should be extended in several directions. First, other 
interest rates should be studied, to know whether the broad features found 
on the U.S. market also holds in other cases. Second, the consequences of 
our modelling for interest rate derivatives should be worked out. Finally, it 
would be interesting to interpret the time dependent smile in option markets 
in term of a 'forward volatility curve', and see how the ideas expressed here 
transpose to this case. We hope to deal with these issues in the near future. 
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